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ABSTRACT 

In order to obtain robust cosmological constraints from Type la supernova (SN la) 
data, we have applied Markov Chain Monte Carlo (MCMC) to SN la lightcurve fitting. 
We develop a method for sampling the resultant probability density distributions 
(pdf) of the SN la lightcuve parameters in the MCMC likelihood analysis to constrain 
cosmological parameters, and validate it using simulated data sets. Applying this 
method to the Joint Lightcurve Analysis (JLA) data set of SNe la, we find that 
sampling the SN la lightcurve parameter pdf’s leads to cosmological parameters closer 
to that of a flat Universe with a cosmological constant, compared to the usual practice 
of using only the best fit values of the SN la lightcurve parameters. Our method will 
be useful in the use of SN la data for precision cosmology. 

Key words: cosmological parameters - cosmology: observations - dark energy - 
supernovae: general 


1 INTRODUCTION 


The use of Type la supernovae (SNe la) as cali¬ 
brated standard candles has led to the discovery of 
dark energy - the acceleratinR expansion o f the Uni¬ 
verse dRiess et al.lll998| : IPerlmutter et allll999ll . The cause 
for the observed cosmic acceleration remains unknown. 
Probing the nature of cosmic acceleration is one of the 
mos t active areas of resear ch today. For recent reviews, 
see Ratra fc Vogelev (|2008ll: Frieman . Turner , fc H u tereil 


200^; Caldwell fc Kamionkowski (l2009l~l: Uzan~ il201 


i tereij 


I Wand ll2010ll : Li et al. ( 2011 1: Weinberg etal.l ~ 2013ll . 

There are major ongoing and planned observational 
projects to illuminate the nature of cosmic acceleration. 
These include the on going Dark Energy Survey (DES) 
dBernstein et al.l l2012lo the dedicated dark energy space 
missi on Euclid, scheduled for launch in 2020 dUaureiis et al.l 
laCTI: and the Large Synoptic Survey Telescope (LSST) 
d Abell et al.ll2009fFI . which is under construction, with first 


light planned for 2019. Dark Energy is one of the main sci¬ 
ence areas for the Wide-Field Infrared Survey Telescope 


* E-mail: mdai@ou.edu; wang@ipac.caltech.edu 
^ http://www.darkenergysurvey.org/ 

^ http://www.euclid-ec.org/ 

® http://www.lsst.org/ 


(WFIRST), w hich could be launched as early as 2023 
dSpergel et aLllioisI l. 

We can expect a dramatic increase in the quantity and 
quality of SN la data in the next decade and beyond. For 
z < 1, thousands of SNe la are expected from DES, and 
hundreds of thousands of SNe la are expected from LSST. 
WFIRST will observe thousands of SNe la at 2 > 1. In 
order to use the SN la data in precision cosmology, it is 
important that we develop robust analysis techniques that 
can be applied to the observed SN la lightcurves. 

In this paper, we apply Markov Chain Monte Carlo 
(MCMC) to SN la lightcurve fitting. We develop a method 
for sampling the resultant probability density distributions 
(pdf) of the SN la lightcuve parameters in the MCMC like¬ 
lihood analysis to constrain cosmological parameters. Both 
of these are new approaches in SN la cosmology. We present 
our methodology in Sec.2, and results in Sec.3. We conclude 
with a summary and discussion in Sec.4. 


2 METHODOLOGY 

As the mechanism of SNe la explosion is still unclear, em¬ 
pirical models are used for fitting SN la lightcurves. The 
lightcurve of an SN la can be characterized by its shape 
and color, which can be corrected to reduce the intrinsic 
dispersion in SN la peak magnitudes. The shape correction 
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utilizes the correlation between SN la peak brightness and 
declin e time (brighter SNe la decline more slowly in bright¬ 
ness) (|Pskovskiil1l977l : iBranchl Il98ll : IPhillipsI 1 19931 . Il999ll . 
The color correction models the variation of SN la colors 
and dust extin ction. _ _ 

Following iBetoule et al.l (201^, we use the SALT2 
m odel (first p ropos ed in IGuv et al.l l^2007^. a nd updated 
in IGuv et ^ ll2010l l and IBetoule et al.l (l2()14h ') for SN la 
lightcurve fitting. SALT2 provides an average spectral se¬ 
quence and color dispersion from the training of a subset 
of SN la lightcurves. For SALT2, The four SN la lightcurve 
parameters are: date of maximum light in the rest-frame 
B band, the amplitude of the spectral sequence (which is 
described conventionally by the peak magnitude in the rest- 
frame B band), lightcurve shape, and color. The date of 
maximum light is a nuisance parameter, since the distance- 
indicator properties of a SN la only depend on its peak 
magnitude, lightcurve shape, and color. The default htting 
procedure is minimizing a to hnd the best ht lightcurve 
parameters from a 4D grid of parameter values. This grid 
method has its limitations: it may result in a local minimum, 
and the error estimate is sensitive to the choice of the grid 
size and the spacing of the grid points. 

In this paper, we use Markov Chain Monte Garlo 
(MGMC) metho d to ht for ligh t curve parameters, using the 
SALT2 model in IBetoule et al.l (l2014l i. The MCMG method 
has the following advantages over a grid-based method: 


(i) The posterior is drawn randomly from the proposed 
distribution; the true probability distribution is recovered 
given enough number of realizations. 

(ii) The statistics can be calculated using multiple chains 
that have converged, leading to robust error estimates. 

(iii) The resultant probability density functions (pdf) are 
smooth and can be utilized in the subsequent cosmological 
analysis. 


2.1 The SALT2 model 

The SALT2 model iIGuv et al.ll2007l : IGuv et al.ll2010l ') is an 
empirical lightcurve model that provides an average spec¬ 
tral sequence and its hrst order variation. The model flux is 
defined as: 

dF 

X) = xqX [Mo{p, X) + xiMi{p, A)-l-...] X exp[cC'I/(A)] 

ClA 

( 1 ) 

Where p is the phase from the maximum light in the rest 
frame, A is the rest-frame wavelength, Mq and Mi are the 
average spectral sequence and its first order variation, GL is 
the color law which is independent of phase. The lightcurve 
parameters xq, xi, and c are determined by lightcurve ht¬ 
ting. To hnd the best ht parameters, the default method 
hnds the minimum using a grid-based method. The is 
expressed as: 

X ~ (^FjYiQflel Fobs) ^SN {.^model Fobs) ^ (2) 

where Fobs is the observed hux after calibration, F^odei 
is the model hux integrated over the observing hlter 
T{X{l + z)) : 

r dF 

Fmodel = il + z) / \—{p,X)T{X{l+z))dX. (3) 


The covariance matrix Csn consists of three parts: a diago¬ 
nal term of the model error, a regular matrix (not diagonal) 
term of color dispersions (K-correction errors), an error ma¬ 
trix (for SNLS) or a diagonal term of the observational hux 
errors (for SN samples other than SNLS): 

Csn ~ F)os,odei C-model Cobs- (4) 

The model parts of the covariance matrix are dependent on 
the lightcurve parameters; they are varied to minimize the 
X^ for a given set of lightcurve parameters. 

We will now discuss in detail the uncertainties in the 
model huxes {Dmodsi term in Eq. (|4]|). The model huxes are 
calculated using Eq. 0 and Eq. 0. The model uncertain¬ 
ties are calculated as 

^model ~ {f^/fiotal) X S X (Vq -f XiVl 2xlVoi) ^ X FjYiodeli 

(5) 

where S, Vo, Vi and Vqi are provided as part of the SALT2 
model, and are dependent on phase and wavelength. Basi¬ 
cally, Vo is the variance in Mq, Vi is the variance in Mi, Vbi 
is the covariance between Mo and Mi. S' is a scale factor, 
/o is the main component of the model hux, which is Mq 
integrated over the observing filter; ftotai is the total model 
hux including the Mi component (not including the color 
term). When ftotai ^ 0, the model uncertainty is set to be 
100 times that of the model hux. Note that in the SALT2 
code, (Vb + xW Vi + 2 * xi * Vbi) is set to 0.0001 when it 
becomes less then zero. 


2.2 MCMC lightcurve fitting 


We use GosmoMC as a generic sampler to generate multiple 
chains using the Metropolis-Bastings algorithm. We assume 
convergence when R — 1 < 0.01 using the G elman and Ru¬ 
bin “R-1” statistic (iBrooks fc Gehn^^ ]j98[l. For a detailed 
discussion on GosmoMC, see iLewis fc Bridle! (l2002l ~). 

To implement SALT2 lightcurve htting using MCMG, 
we hrst need to understand how the SALT2 default grid- 
based method works. In particular, how the covariance ma¬ 
trix from Eq. 0 is handled, since the contribution from 
model uncertainties depends on the values of the lightcurve 
parameters (which are being htted). There is an optional 
“update weights” feature in the p ublic SALT2 code, but 
it seems not t o be used by either iGonlev et al.l (1201 ill or 
IBetoule et al.l ll2014ll : we are only able to reproduce their 
results without using this option (see SeclSj. However, by 
not updating the weights for the covariance matrix in the 
SALT2 code, it does not mean keeping the covariance ma¬ 
trix hxed. It means deriving a converged covariance matrix 
by doing the following: 


(i) Initial ht with xi Hxed to be 0, without including 
model uncertainty. 

(ii) Second ht with xi allowed to vary, without including 
model uncertainty. 

(iii) Iterations of fits with model uncertainty included, 
until the changes in all parameters are less than 0.1 times 
the errors in the parameters. 


At the beginning of each iteration in step 3, the co- 
variance matrix is recalculated using the parameters from 
the previous step or iteration. And the covariance matrix 
is kept fixed during this iteration of the fit. In MCMG we 
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cannot do the fit by steps as in the grid-based method, so 
we calculate the covariance matrix using fiducial lightcurve 
parameter values and keep it unchanged during the MCMC 
lighcurve fitting. We assume that likelihood oc exp(—x^/2), 
with defined in Eq. ©, and carry out the following steps: 

(i) Perform a grid-based fit to obtain a set of fiducial val¬ 
ues of lightcurve parameters (which are the best fit values 
from this fit). 

(ii) The covariance matrix Csn is calculated using the 
fiducial values of the model and lightcurve parameters. 

(iii) An MCMC likelihood analysis is performed to ob¬ 
tain the lightcurve parameters, while fixing the covariance 
matrix to be that calculated using the fiducial lightcurve 
parameters. 

A compelling reason for us to fix the covariance matrix 
for MCMC lightcurve fitting to be that calculated at fiducial 
values of the lightcurve parameters is as follows. The model 
uncertainty from Eq. is dependent on the lightcurve pa¬ 
rameters (especially xi and time of maximum flux), so it is 
varied as the lightcurve parameters are vared. For example, 
if xi is a very large number, it is possible that ftotai will 
become negative, which means that the model uncertainty 
is set to 100 times that of the model value, leading to a very 
small — a numerical artifact that can bias the lightcurve 
fitting. We have found that simply updating the covariance 
matrix as the lightcurve parameters vary during the MCMC 
steps can lead to unreasonable results for the lightcurve pa¬ 
rameters, since the model errors can be much larger when 
the lightcurve parameters fall out of the reasonable ranges. 

To assess the implications of our choice of fixing the co- 
variance matrix for MCMC lightcurve fitting, we have com¬ 
pared the lightcurve parameter results using two different 
fiducial values — one set contains the best fit results from a 
grid-based fit, the other set the best fit results from a grid- 
based fit without model errors. The lightcurve parameters 
we get from MCMC using the two different covariance ma¬ 
trices are very similar to each other (with a mean difference 

10““^ for all three parameters used in a cosmological fit, 
i.e. ms, xi, and c), and lead to almost identical cosmologi¬ 
cal constraints. This is not surprising, since the cosmological 
constraints are not very sensitive to the lightcurve param¬ 
eters, as long as no erroneous lighthcurve parameter values 
are used. We conclude that it is reasonable to fix the covari¬ 
ance matrix at fiducial values of the lightcurve parameters 
in MCMC lightcurve fitting. 


2.3 Cosmological analysis 

Having derived SN la lightcurve parameters using MCMC, 
we can use them to derive cosmological constr aints. We fol¬ 
low th e definition for the model magnitude in IConlev et al.l 

(l201lli : 

uimod = SlogjQ - oxi-I-/3 c-I-Af, (6) 

where Dl is a redefined luminosity distance that is inde¬ 
pendent of the Hubble constant, a, /3 and M are nuisance 
parameters which describe the shape and color corrections 
of the lightcurve, and the SN absolute magnitude in com¬ 
bination with the Hubble constant. In order to model the 
dependence of SN la intrinsic brightness on the host galaxy 


mass, M is defined as a function of the host galaxy stellar 
mass Mhost (in units of solar masses): 


Ad = / Mhost < 10 

\ M 2 for logiQ Mho3t > 10 

The x^ is then 

X^ = Am^ • ■ Am, 

where Am — niB — mm od, and mg is calcu lated by ms = 
—2.5 logjQ(a;o) + 10.635 (iMosher et al.l[20l3 ~l. C is the co- 
variance matrix. We use the same co variance matrix in ou r 
cosmological analysis as that used bv iBetoule et al.1 ll2014h . 
Note that rumod is defined in Eq.®. 

In this paper we assume a fiat Universe with constant 
dark energy equation of state w, since SN la data alone do 
not provide meaningful constraints on additional cosmolog¬ 
ical parameters. The Hubble constant free luminosity dis¬ 
tance Hi is defined as: 

Hi =c-'Ho(l + «hei)r(2), (9) 

where 2hei is the heliocentric redshift, 2 is the CMB-frame 
redshift (i.e., the cosmological redshift), r(z) is the comoving 
distance: 

riz) = cHo^r{z) (10) 


(7) 

( 8 ) 



( 11 ) 


E{z) = H{z)/Ho (12) 

With the assumption of flat Universe and constant w, 

H{z) = -= Hox/n^(l -b 2)3 -b Odb)! -b 2)3(1+-) ( 13 ) 

a V 

and flm + Qde = 1 (the radiation contribution is negligible 
here). 


2.4 pdf sampling 


The MCMC analysis gives the marginalized pdf’s of the 
lightcurve parameters. Those pdf’s contain the distribution 
and error information of the lightcurve parameters and can 
be used in a cosmological analysis. We have developed a 
method to sample the pdf’s, and derive cosmological results 
by combining results from different sets of lightcurve param¬ 
eters drawn from the pdf’s, as described below. 

For each SN, we choose N points with equal probabil¬ 
ity intervals from the pdf of each lightcurve parameter, xo, 
xi and c, with probabilities equal to Pi, P 2 , ... , Pjv- This 
gives sets of lightcurve parameters, with each set having 
the same SNe with different values of lightcurve parameters. 
We then use these sets of lightcurve parameters to fit cos¬ 
mology, resulting in sets of cosmological parameters. We 
combine the results in the following way to get the combined 
cosmological parameters: 


PiPjPk Sijk 

'^i,j,k PiPjPk 


(14) 


where s are the cosmological parameters from sampling the 
pdf’s of the lightcurve parameters, Sijk are the cosmolog¬ 
ical parameters derived from the data set with a given 
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set of lightcurve parameters drawn from the pdf’s of the 
lightcurve parameters. The cosmological parameters from 
different data sets are weighted by the product of the rela¬ 
tive probabilities of the three lightcurve parameters, Pi, Pj, 
and Pk- Using only lightcurve parameters that correspond 
to the peaks of the pdf’s gives Pi = Pj = Pk = 1; this is sim¬ 
ilar to the usual practice of using only the best fit lightcurve 
parameters in the grid-based method for SALT2 (which is 
equivalent to using the peaks of the mean likelihood distri¬ 
butions from MCMC). 

It is not practical to densely sample the pdf’s of the 
lightcurve parameters. In order to gauge how the cosmolog¬ 
ical parameter constraints depend on the sampling density 
of the lightcurve parameter pdf’s, we study the following 
cases: 

(i) N=3. We choose three points from the pdf, with the 
probabilities Pi = 1 and P 2 = P 3 = 1/2. This results in 3 
points on the pdf: the peak, and the half height point on 
either side of the peak. 

(ii) N=7. We choose seven points on the pdf. Pi = 1 (the 
peak), P 2 = P 3 = 1 / 2 , P 4 = P 5 = 3/4, and Pe = P 7 = 1/4. 
This divides the pdf in 1/4 segments in height, resulting in 
three points on either side of the peak. 

(iii) N=15. In this case 15 points are chosen from the 
pdf, the probabilities are 1, 1/2, 3/4, 1/4, 7/8, 5/8, 3/8, 
1 / 8 . 

(iv) N=19. When dividing the pdf into 10 equal proba¬ 
bility intervals with probabilities Pi = i/10{i = 1, 2,..., 10), 
we get 19 points on each pdf. 

We will show that pdf sampling is converged with increased 
N in Section EH 


3 RESULTS 

We have applied our methodology for SN la lightcurve fit¬ 
ting using MCMC and constraining cosmology with sam¬ 
pling the pdf’s of SN la lightcurve parameters to the 
“Joint Lightcurve A nalysis” (JLA) data set of SNe la from 
iBetonle et al.l (l2014ll . which combines the SNLS and SDSS 
data of SNe la in a consistent, well-calibrated manner. We 
do our MCMC lightcur ve fitting using the c alibrated photo¬ 
metric data provided bv iBetoule et al.l (l2014h . The JLA sam¬ 
ple, a s an extension to the Cll Compilation dConlev et al.l 
Im3), contains a combination of data sets of 740 spec¬ 
troscopically conf irmed SNe la from several low-z samples 


(2 < 0.1) (mostlv Hicken et al.l 20091: Contreras et al. 

l20inl: 

iHamuv et al.lll996l: iJha et al.ll2006l: iRiess et al.lll999 

), the 


full three-year SD SS-II supernova survey (0.05 < 2 < 0.4) 
dSako et ai.ll20lj ’). the hrst thre e year s data of the SNLS 
survey (0.2 < 2 < 1) lICuv et al.ll20ld : IConlev et al.ll201ll ~l 
and a coupl e of high redshift HST SNe (0.7 < 2 < 1.4) 
dRiess et al.ll2007t l. The photometry of SDSS and SNLS is 
recalibrated. The SALT2 model is retrained using the joint 
data set. 

3.1 Definitions 

We will show comparative results on cosmological con¬ 
straints using different sets of lightcurve parameters we have 
obtained in different approaches. These are: 


(i) S ALT2 : This is the p ublished lightcurve parameter 
set from iBetoule et ^ |201J), we use it directly in the cos¬ 
mological analysis as the base of comparison to other sets. 

(ii) GRID-SALT2 : We obtain this set of lightcurve pa¬ 
rameters by running the published version of the SALT2 
CO defl adding in the bias correction term to the peak mag¬ 
nitude and the uncertainties in redshift, lensing and intrinsic 
dispersion to the magnitude uncertainties. All the values of 
the terms above are the same as tho se used in the JLA set 
as described in IBetoule et al.l d2014ll . This is supposed to 
reproduce the results of SALT2. All other sets of lightcurve 
parameters described in the following are processed as de¬ 
scribed here. 

(iii) GRID : We use our own grid-based code using the 
SALT2 model to calculate the this is an important cross¬ 
check, to ensure that we understand all the nuances of the 
public SALT2 code and its output. We will use this as the 
grid method to compare with our MCMC analysis. We ob¬ 
tain this GRID set of lightcurve parameters by minimizing 
the using exactly the same approach as the SALT2 code 
does (he. using the function minimization package called MI- 
NUITlj) . We expect to get the same values for the lightcurve 
parameters within numerical errors. 

(iv) MCMC-LIKE : The MCMC chains can be used to 
calculate mean likelihood distributions; its maximum corre¬ 
sponds to the least value from the grid method. If the 
grid method ever falls into a local minimum, the maximum 
likelihood value of the MCMC chains would instead give the 
correct global minimum upon convergence. 

(v) MCMC-MARGE : The MCMG chains can be used 
to obtain marginalized one-dimensional distributions of the 
fitted parameters; these give the standard error distribu¬ 
tion information from an MCMG analysis. For this MCMG- 
MARGE set we use the means of the marginalized pdf’s 
as lightcurve parameters for htting cosmology (with no pdf 
sampling). We will discuss the pdf sampling results in Sec¬ 
tion EH 

In an MCMG analysis, the differences of the marginal¬ 
ized and the mean likelihood distributions indicate non- 
Gaussianity, although it is possible to have a non-Gaussian 
distri bution where both curves are the same llLewis fc Bridlel 
I 2 OO 2 II . In general, the marginalized pdf differs more from 
the mean likelihood for parameters that are less well con¬ 
strained by the data. Mean likelihood shows how good a 
ht you could expect if you drew a random sample from the 
marginalized distribution. It is customary to quote marginal¬ 
ized constraints in a cosmological analysis. We will follow 
this practice in this paper regarding the lightcurve parame¬ 
ters, except when we need to make a direct comparison with 
the grid-based method (which gives results that are equiv¬ 
alent to the mean likelihood). However, when showing the 
cosmological constraints, we give the mean likelihood pdf’s 
instead, as our tests with simulated data sets show that the 
peaks of the mean likelihood pdf’s of the cosmological pa¬ 
rameters are less biased than the marginalized means. (For 
more details, see section 13.4.Il l 


http://supernovae.in2p3.fr/salt/doku.php 
® http://seal.web.cern.ch/seal/snapshot/work- 
packages / mathlibs / minuit 
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Figure 1. Mean likelihood pdf’s of the parameters for the cos¬ 
mological fit for the grid sets, using all 740 SNe from the JLA 
data set. Black solid lines represent the GRID set; red dashed 
lines represent the GRID-SALT2 set; blue dotted lines are for the 
SALT2 set. 

3.2 Reproducing the SALT2 Results 

We first compare results of the three grid sets - SALT2, 
GRID-SALT2 and GRID, to verify that our grid-based code 
is correct. In this cross-check exercise, we found that we 
are only able to reproduce the lSetoule et ahl (l2014h results 
without using the “update weights” option in the public 
SALT2 code, but instead carrying out the steps as described 
in Sec l2.2l 

The lightcurve parameters from the three sets are very 
similar. Fig[T] shows a comparison of the cosmological con¬ 
straints from SALT2, GRID-SALT2, and GRID. As ex¬ 
pected, all three approaches give nearly the same constraints 
on the cosmological and SN nuisance parameters. 

3.3 MCMC vs GRID Comparison 

We now compare the results of MCMC and GRID. We find 
that the lightcurve parameters from the MCMC-LIKE set 
are also consistent with that from the GRID set. However, 
the lightcurve parameters of the MCMC-MARGE set are 
offset from the MCMC-LIKE (or the GRID) set for some 
SNe. The cosmological results are shown in Fig. [2] The 
MCMC-LIKE set gives similar cosmological constraint as 
the GRID set, as expected. The MCMC-MARGE set shows 
pdf’s shifted from the MCMC-LIKE (or the GRID) set in 
Q.m and in; this is not surprising given the difference between 
the marginalized and the mean likelihood in an MCMC anal¬ 
ysis (see the discussion near the beginning of Sec[3]). The 
nuisance parameters are well constrained and have similar 
constraints in all sets. 

The peaks of the mean likelihood pdf’s and their 68% 
confident intervals of the parameters are listed in Table [T] 
We find the approximate 68% confidence levels by finding 
the parameter values where the probability has dropped by a 
factor of Note that for some sets there is only one up- 


Figure 2. Mean likelihood pdf’s of the parameters from the cos¬ 
mological fit for the MCMC sets comparing with the GRID set, 
using all 740 SNe from the JLA data set. Black solid lines rep¬ 
resent the MCMC-MARGE set; red dashed lines represent the 
MCMC-LIKE set; blue dotted lines represent the GRID set. 


per limit in Q.m, because the pdf’s are truncated at firn = 0 
before the probability has dropped by a factor of The 

grid methods (SALT2, GRID-SALT2, GRID) show similar 
values and errors. The MCMC-LIKE has consistent values 
with the grid methods but has slightly higher errors in flm 
and w. The MCMC-MARGE set show even higher errors in 
w, but firn and w are closer to the “concordance model” of 
fim = 0.27 and w = —1. 


3.4 Cosmological Constraints from Sampling the 
pdf’s of SN la Lightcnrve Parameters 

We now present cosmological constraints derived from sam¬ 
pling the pdf’s of SN la lightcurve parameters. For this work, 
we have to limit our analysis to SNe la with well-behaved 
lightcurve parameter pdf’s. For most SNe, the pdf’s of 
their lightcurve parameters are well-behaved, single-peaked 
smooth bell curves. However, there are several exceptions 
with multi-peak pdf profiles for the lightcurve parameters. 
We have tracked the multi-peak profiles to data quality is¬ 
sues: some SNe have no data after the maximum light, some 
SNe have lightcurves with too few data points, or very noisy 
data. We exclude those problematic, multi-peak SNe from 
our cosmological analysis. We also exclude other SNe that 
don’t have any data in any bandpass after the peak magni¬ 
tude as a quality cut. This results in a set of 729 SNe la; we 
will use only this set in our analysis from this point on. 

The effects of excluding the 11 problematic SNe (listed 
in Table[2]) from the cosmological analysis are shown in Fig. 
[3] The constraints on fim and w are noticeably shifted by 
excluding these 11 SNe, indicating that cosmological results 
could be biased by including poor quality data. We observe 
the same effect both using the original JLA data and using 
our MGMC-fitted lightcurve parameters. 
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Table 1. Parameters from the cosmology fit using all 740 JLA SNe 




w 

a 

/3 

Ml 

M2 

SALT2 

0.196+°“® 

-0.795l°i®i 

n 141+0008 

Q -| r\'7+0.078 

O.iu/ -0.088 

24.118l°;“I 

Z4i.UOU_0.029 

GRID-SALT2 

0.185+°+°° 

-0.806t°J®8 


o.J.u»_o 090 

24.1181°°°® 

Z4i.u4io_o 028 

GRID 

0 i 7 i +°.145 

- 0 . 769 t°oill 

n 141+0006 

Q 1 r\ir+0.085 
O.iUO_0 070 

24.1171°;°°^ 

Z4i.U4i4t_o 025 

MGMG-LIKE 

0.171+°.157 

rv <7 qq~I" 0-257 

U. <oo_ 0.282 

n 141+0006 

'J--L^-L_0.007 

Q nss+0'087 

CS.UOO_0.082 

z^.izo_o.o32 

Z4i.uoU_0.030 

MCMC-MARGEt 

U.Z40_q jgj 

-0.87lt°;2°| 

0 141+0.007 
'J--L^-L_0.007 

o.j-uu-o 083 

24 110+0025 

^‘i.Tiy-O 029 

24 045+0 029 
Z4i.u4io_o 030 


t The MCMC-MARGE set uses the marginalized means of the lightcurve parameter pdf’s. 


Table 2. Problematic SNe 


SN name 

1 

Lancaster 

2 

Patuxent 

3 

SDSS 11206 

4 

SDSS14318 

5 

SDSS16619 

6 

SDSS16737 

7 

SDSS16793 

8 

SDSS19067 

9 

SDSS21510 

10 

SDSS21669 

11 

Torngasek 



0 0.2 0.4 



0.12 0.13 0.14 0.15 0.16 

cx 



24.05 24.1 24.15 24.2 

Ml 



24 24.05 24.1 

M2 


I-729SN-740Sn| 


cosmological model, added with and ran¬ 

domly drawn from the lightcnrve parameter pdf’s we obtain 
using our MCMC lightcurve fitter, and a Gaussian scatter 
(p = 0,(7^ = 0.1^): 

(15) 

The nuisance parameters a, /3 and M are defined in Eq. |6l 
and are fixed with fiducial values. 

We generate 1000 sets of simulated data and perform 
cosmological analysis to each set of data, with and without 
pdf sampling. When applying pdf sampling to the simulated 
dataset, we only sample the xi and c parameter, and only 
3 points are chosen from each pdf to speed up the process. 
We list the input parameters and the means and standard 
deviations of the 1000 sets of resultant cosmological param¬ 
eters in Table [3] We have shown both the peak values of 
the mean likelihood pdf and the means of the marginalized 
pdf. The two values and their standard deviations show dif¬ 
ferences in the Q.m and w parameter, which is due to the 
differences in the meaning of the two kind of pdf’s. We have 
briefly discussed the differences between the mean likelihood 
pdf and the ma rginalized pdf in Sec. 13.11 For more detailed 
discussions, see iLewis fc Bridle! (l2002ll . Since the peaks of 
the mean likelihood pdf are less biased in general, we only 
show the mean likelihood pdf’s for cosmological constraints 
in this paper. When quoting the peak values of the mean 
likelihood pdf’s, comparing to not using pdf sampling, ap¬ 
plying pdf sampling gives Q.m and w that are closer to the 
input parameters, which is expected as pdf sampling uti¬ 
lizes more information from the lightcurve parameter pdf’s. 
However, the standard deviations are larger in these two pa¬ 
rameters when applying pdf sampling, as the values spread 
in larger ranges. 


Figure 3. Mean likelihood pdf’s of the parameters from the cos¬ 
mological fit for different numbers of SNe la using the MCMC- 
MARGE set. Black solid lines show the results from 729 SNe 
excluding the 11 problematic SNe; red dashed lines represent the 
whole SNe sample including all 740 SNe from the JLA data set. 

3.4.1 Tests with simulations 

Before showing the results of pdf sampling applied to the 
JLA dataset, we first show results from the simulated 
datasets to illustrate that pdf sampling gives less biased cos¬ 
mological results than the usual practice without pdf sam¬ 
pling. To generate the simulated datasets, we replace the 
bias-corrected B-band peak magnitude (ms) in the JLA 
dataset with true peak magnitude calculated from a fiducial 


3.4.2 pdf sampling with the JLA dataset 

We now proceed to implement pdf sampling of the lightcurve 
parameters in our cosmological analysis, as described in 
Sec l2.4l to the set of 729 SNe la, as described below: 

(i) PDF-COMBINED-3 : We draw 3® = 27 sets of 
lightcurve parameters from the pdf’s as described in Sec¬ 
tion We also apply the same bias correction as used in 
the JLA set to each individual set of lightcurve parameters. 
The cosmological results of the individual sets are combined 
as described in Section[2]4]to obtain the cosmological results 
with pdf sampling of lightcurve parameters. 

(ii) PDF-COMBINED-7 : Similarly, 7^ = 343 sets of 
lightcurve parameters are drawn from the pdf’s and the cos¬ 
mological results are combined. 
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Table 3. Test with simulated data 




w 

a 

/3 

Ml 

M2 

input values 

0.3 

-1.0 

0.14 

3.1 

24.11 

24.04 

w/o pdf (like®) 

0.205±0.098 

-0.928±0.196 

0.152±0.006 

3.543±0.088 

24.114±0.019 

24.034±0.021 

w/ pdf (like^) 

0.258±0.111 

-0.990±0.288 

0.152±0.007 

3.541±0.089 

24.123±0.032 

24.038±0.036 

w/o pdf (marge^) 

0.247±0.057 

-1.038±0.163 

0.152±0.006 

3.548±0.088 

24.112±0.018 

24.032±0.020 

w/ pdf (marge^ 

0.287±0.055 

-1.121±0.181 

0.148±0.006 

3.545±0.089 

24.111±0.018 

24.029±0.020 


^ The means and standard deviations of cosmological parameters of the 1000 simulated datasets, by quoting the peak 
values of the mean likelihood pdf’s of the cosmological parameters; 

^ The means and standard deviations of cosmological parameters of the 1000 simulated datasets, by quoting the means 
of the marginalized pdf’s of the cosmological parameters. 


(iii) PDF-COMBINED-15 & PDF-COMBINED-19: 
15® = 3375 and 19® = 6859 sets of lightcurve parameters 
are drawn respectively. These ensure that the combined 
results have converged with increasing N. (See Fig. ^ 

We compare the results from using pdf sampling of 
lightcurve parameters with the GRID set and MCMC- 
MARGE (no lightcurve parameter pdf sampling), shown in 
Fig. [5] Since we have already shown that the pdf sampling 
results are converged with increasing N, we only show one 
set of results from using pdf sampling - PDF-COMBINED- 
15. 

Note that the pdf’s for firn and w are shifted, compared 
to the results using the GRID set (grid-based method with 
no pdf sampling). When comparing with the results using 
only the marginalized mean values of the lightcurve param¬ 
eters (without pdf sampling), the pdf’s are also shifted, and 
pdf-sampling gives a little tighter constraints. The peak of 
the mean likelihood pdf’s and their 68% confidence inter¬ 
vals of the cosmological analysis using 729 SNe excluding 
the problematic ones are shown in Table [4] When excluding 
the problematic SNe, we get slightly larger values, and 
the w values are closer to —1. In Fig. [SI we show the corre¬ 
sponding 2D mean likelihood contours of the fitted parame¬ 
ters of PDF-GOMBINED-15, GRID and MGMC-MARGE. 
It is interesting to note that pdf sampling of SN la lightcurve 
parameters, and using the means of the lightcurve param¬ 
eter pdf’s, leads to cosmological constraints closer to a flat 
Universe with a cosmological constant. 



Figure 4. Mean likelihood pdf’s of the parameters from the cos¬ 
mological fit with sampling of the SN la lightcurve parameter 
pdf’s, using 729 SNe la from the JLA data set (excluding 11 
problematic ones). Black solid lines are results from sampling 3 
points on each pdf (PDF-COMBINED-3); red dashed lines are 
results from sampling 7 points on each pdf (PDF-COMBINED- 
7); blue dotted lines are results from sampling 15 points on each 
pdf (PDF-COMBINED-15); magenta dash-dotted lines are re¬ 
sults from sampling 19 points on each pdf (PDF-COMBINED- 
19). 


4 SUMMARY AND DISCUSSION 

We have developed a method to utilize the probability den¬ 
sity distribution (pdf) of SN la lightcurve parameters in 
cosmological analysis using SNe la. First, we have applied 
Markov Ghain Monte Garlo (MCMG) to SN la lightcurve 
fitting, in order to obtain smooth and well-behaved pdf’s of 
SN la lightcurve parameters. Then we derived cosmological 
constraints with sampling of the pdf’s of the SN la lightcurve 
parameters. For a complementary approach of sampling the 
un derlying SN la popu lation in cosmological model fitting, 
see lMa,rch et al.l (l201ll) . 

In order to validate our method, we applied it to 1000 
sets of simulated SN la data. We found that compared to not 
using pdf sampling, applying pdf sampling gives Qm and w 
that are closer to the input parameters, which is expected as 
pdf sampling utilizes more information from the lightcurve 


parameter pdf’s. We also found that the peak values of the 
mean likelihood pdf’s are closer to the input values than the 
marginalized means. 

Our method differs from the usual approach in two 
ways: (1) We use MGMC, instead of a grid-based method, in 
fitting the SN la lightcurve parameters; (2) We sampled the 
pdf’s of the SN la lightcurve parameters in the cosmological 
analysis, instead of just using the peaks of the pdf’s. 

We have applied our method to the Joi nt Lightcurve 
Analy sis (JLA) data set of SNe la derived bv iBetoule et ah] 
ll20l4 l. which combines the SNe la from SDSS and SNLS in 
a consistent, well-calibrated manner. Interestingly, we find 
that the resultant cosmological constraints are closer to that 
of a flat Universe with a cosmological constant, compared to 
the usual practice of using only the best fit values of the SN 
la lightcurve parameters. 

The JLA set has a bias correction term which is deter- 
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Table 4. Parameters from the cosmology fit using 729 JLA SNe (excluding 11 problematic ones) 




w 

a 

p 

Ml 

M2 

SALT2 

n 

U.ZOO_Q JgQ 

-0.848t°;f7° 

U.i4±U_o 006 

o.uy^_Q Q<^2 

^^.Tzo_o 024 

24 054+0-028 

Z4i.uO4_0 030 

GRID-SALT2 

0.258l“;i°4 

-0.866t°i®° 

0 141+0-007 

^•-L‘±-L_0.006 

o.uyu_Q Qg3 

24 123+0-026 
^‘i.TZO_o 028 

Z4i.u01_o 030 

GRID 

0.249+°“® 

-0.830t°nlm 

0 141+0-007 

'J--L^-L_0.006 

3 093+0-083 

o.uyo_Q Qg7 

24 124+0-024 

24 053+0-026 
Z4l.UOO_o 031 

MGMG-LIKE 

u.zuo_Q 205 

-0-85+0.359 

0 141+0-007 
^--L‘±-L_o.o07 

q nyi +0.083 
o.ui -*-_o 080 

24.1261°°®° 

24 049'*'^'^^® 
Z4i.u^y_0 031 

MCMC-MARGEt 

U.OTZ_q jgg 

-0.965l°®°° 

0 142+0-006 
u.T^z_o 007 

o.uyu_o og7 

24.1281°°“ 

24 040"*”^'^^^ 
Z4l.U^U_o 024 

PDF-COMBINED-3 

0 339+0 065 

u.ooy_Q jgj 

-0-910io.339 

0 142+0-006 
u.T^z_0 005 

3.0761°°®® 

24 124+0-023 

24 046'*'^'^^^ 
Z4i.u^O_o 039 

PDF-COMBINED-7 

n 99 c:;+ 0-072 
U.OCiO_Q 

-0.975l°+“ 

0 142+0-005 

u.T^z_o 006 

3.0981°°°° 

24 107+0-031 
Z41.1U/ _o 026 

Z4l.UZO_o 023 

PDF-GOMBINED-15 

n Qoci+O-OSi 

U.OZO_o ig2 

U.yiO_Q 344 

0 143+0 005 

U.T^O_0 QQY 

3.0841°;°^® 

24 110+0-026 
Z4i.iio_0 032 

24 036'*'^'^^^ 
Z4i.uOO_o 028 

PDF-GOMBINED-19 

u.ouy_Q jg2 

—0 912+0-343 

u.yi^_Q 330 

0 143+0 006 
U.T^0_Q QQg 

3.0961°°^° 

24 110+0-028 

Z4l.liO_o 029 

24 035+0-034 

Z4i.uoo_o 026 


f The MCMC-MARGE set uses the marginalized means of the lightcurve parameter pdf’s. 



Figure 6. Joint confidence level contour plots from the cosmological fit with and without sampling of the SN la lightcurve parameter 
pdf’s, using 729 SNe la from the JLA data set (excluding 11 problematic ones). The contours are 68% and 95% confidence levels. Thick 
black solid contours are are results from sampling 15 points on each pdf (PDF-COMBINED-15); thin red solid contours are results from 
using the GRID set; blue dotted contours are results from using the MCMC-MARGE set, without pdf sampling. 


mined using simulat ions generated by the SNANA software 
llKessler et al.l l2009l ) , by comparing the reconstructed dis¬ 
tance using SA LT2-fitted lightcurve parameters to the sim¬ 
ulation inputs (iBetoule et al.ll2014^ . We use this bias cor¬ 
rection in our analysis, assuming that it is independent of 
the fitting technique. Ideally we should perform our own 
bias calculation by fitting the simulations with our MCMC 
lightcurve fitter. We will leave this for future work. We also 
note that the significantly smaller marginalized errors (com¬ 
pared to likelihood errors) from the simulated datasets could 
be due to the simulated datasets not being sufficiently real¬ 
istic. We will investigate this further in future work. 

As SN la data increases in both quantity and quality, 
our method will be useful in the quest to illuminate the 
nature of dark energy using cosmological data. 
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